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QUANTUM PHASE TRANSITION FROM A STRANGE METAL 

TO A FERMI LIQUID IN THE NORMAL STATE OF 

CUPRATE SUPERCONDUCTORS 

ANDREW DAS ARULSAMY 



Abstract. The normal state metallic behavior above the superconducting 
fl\ transition temperature (Tsc) remains one of the least understood properties 

r^) since its discovery in cuprates by Bednorz and Miiller [J. Apart from the 

pseudogap (T*) phenomenon above Tgc, there is this T-linear resistivity that 
OO was first recognized by Anderson [2] as strange simply because it does not 

Cm obey the low-temperature T^ electron-electron scattering rate predicted from 

the Fermi liquid theory. Here, we prove that the finite temperature continu- 

' I ' ous quantum phase transition (jSQPT) is responsible for the transition from a 

^ strange metal to a Fermi liquid for a given chemical composition. Our proof 

' follows the Green function formalism and the energy-level spacing renormal- 

•4.^ ization group method. Interestingly, the ^QPT stated above can be related 

^ to the proof derived by Parameswaran, Shankar and Sondhi [3] , in which, the 



one-to-one correspondence between the Fermi gas and Fermi liquid is indeed a 
special case. The theory is further exploited to give unambiguous explanations 
r^ as to why and how doping gives rise to spectral weight transfer observed in the 

I soft X-ray absorption spectroscopic measurements for the La-Sr-Cu-O system. 

o 

O PACS: 74.72.-h, 71.10.Hf, 74.25.F-, 74.25.Jb 



1. Introduction 



Landau's Fermi liquid theory (FLT) was originally developed to explain the liquid 
state of Helium-S. Later, it was applied to explain the free-electron metals (with 
1/-^ one-to-one correspondence between FLT and free-electrons), where it has fulfilled 

^^ its purpose with unprecedented accuracies such that its predictions are qualitatively 

r^ comparable with experiments [H |SJ [SI [71 |H] ■ However, the knowledge derived from 

f^ FLT is found to be inadequate (not inaccurate) to explain certain metals with 

'~H strongly interacting electrons. For example, the metallic property of heavy fermions 

'^ j and high temperature superconductors above the critical transition temperature do 

J> not obey FLT [HI [lUj- When one comes to think of it, this is not surprising at all 

because the metallic behavior of strongly-interacting electrons, physically should 
not obey the non- interacting fermions (or Fermi gas). On the other hand, any 
C^ new theory put forward to explain the metallic property in the presence of strong 

interaction(s) will have to be inadequate to explain the free-electron metals, which 
is a necessary condition. This condition is necessary because removing or adding 
this interaction Hamiltonian {Hint{x)) is always by hand, even if it is done with 
appropriate limits {Hint{x) — >■ or x — > oo), where x is some parameter inversely 
proportional to the interaction strength. Such theoretical inadequacies (between 
FLT and the prospective theory) imply the existence of a quantum phase transition 
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due to changing chemical composition via doping (above and/or below a certain 
critical temperature). Details on the zero temperature quantum phase transition 
(okQPT) due to changing pressure, doping and magnetic fields, and Mott-Hubbard 
transition are given in Refs. [H HI HI HI [HI HI [JT] . 

There are other types of critical temperatures such as the superconducting (and 
the pseudogap) transition temperatures (Tgc and T*) for cuprates, the Curie tem- 
perature (Tc) for ferromagnets, the Neel temperature (Tn) for antiferromagnets, 
or the ferroelectric transition temperature (Tfc) for ferroelectrics, in addition to 
the metal-to-insulator transition temperature. Here, a new critical phenomenon 
is proven to exist, which can be associated to the continuous finite temperature 
quantum phase transition (^QPT). For example, from a strongly interacting metal 
(a strange metal) to a Fermi metal. However, there is an important exception, in 
which, the ^QPT developed and proven in this work does not include the super- 
conductor transition temperature and the pseudogap phase in any way. As such, 
the 2^QPT must be above (not below) the T* if Tsc < T* . The term "non-metals" 
means semiconductors or insulators such that the electrons in the conduction or 
overlapped band cannot be treated as free- electrons. The superconducting phase 
and the regime in its vicinity above Tgc has been deliberately excluded because 
we are not sure (beyond any doubt) how the energy level spacing of a given com- 
pound evolve with chemical doping to give rise to superconductivity (this work is 
in progress). 

In the early days of the ionization energy theory (lET), it is true that we have 
developed resistivity models for cuprates for all temperatures between Tgc (includ- 
ing T*) and room temperature using exotic particles (due to charge-spin separa- 
tion [31 HSl HH HO]) J namely, spinous (uncharged electrons) and holons (spinless 
charged bosons) [5T[ [52]. A proof is also developed (using the loffe-Larkin con- 
ductivity rule j23j) to justify the above models by invoking the co-existence of 
electrons, spinous and holons |24| . However, the theory developed and proven here 
is entirely based on the strongly correlated electrons within lET and the Green 
function formalism, without requiring such exotic particles. For the records, lET 
is not a non-Fermi liquid theory. In particular, these strongly correlated electrons 
require certain additional conditions compared to the well-known Fermi liquid the- 
ory and the Fermi gas concept. In view of this, indeed, we have been unfaithful and 
have switched side (back to electrons) after the discovery of (i) the pseudogap-like 
features in manganites |2S] and (ii) the existence of ^ 7^ even in the presence of 
crossed energy levels [26j . To make sure the readers and I are in resonance, a new 
Section (Further analysis) is added toward the end, to iron out any unclear issue 
that may have left unexplained. 

Our main conclusions are as follows; (a) the phase transition from a strange 
metal to the Fermi liquid is due to the finite-temperature continuous quantum 
phase transition (QPT) as a result of the insignificant effect from the energy-level 
spacing, (b) Fermi gas is a special case having a one-to-one correspondence with 
the Fermi-liquid if the energy-level spacing is zero, (c) the spectral weight trans- 
fer due to doping in cuprates observed in the soft X-ray absorption spectroscopic 
measurements can be unequivocally captured with the Green function developed 
within lET, and (d) the emergent phenomenon put forth by Anderson |27| can be 
related to quantum phase transitions. 
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2. Renormalized Green function 

Before we explore the electronic phase transitions in superconducting cuprates 
above Tgc, we first need to renormalize the relevant interaction term in the Hub- 
bard Hamiltonian via the energy-level spacing renormalization technique within the 
Green function formalism. The reason for doing this is to expose any interesting 
electronic phase transition due to crossed energy-levels, giving rise to overlapped 
bands such that ^ 7^ 0. However, the free-electron metallic phase is well-known to 
exist in the overlapped region, which satisfies the weakly-interacting Fermi liquid 
that has a one-to-one correspondence with the Fermi gas when the Mott-Hubbard 
gap (Uh) closes and ^ = 0. The conditions Uh = and ^ 7^ can be understood by 
noticing that the energy-level crossings at a certain k-point, say at ki may imply 
EaCki) — Ei,(ki) but ^ 7^ due to interactions where Ea and Ei, are the eigenvalues 
for the solved two-level Hamiltonian 

il(k)(^,(k) = [ifo(k) + V{k)](^,(k) = [/i,(k) +V,(k)]^a(k) = Ea{k)^a{k), 

(2.1) i/(k)(^,(k) = [Ho{k) + V{k)]M^) = [hb{k) + vb(k)]^fc(k) = Eb{k)M^). 

Here, i?o(k) is the basic Hamiltonian excluding all interactions, whereas V(k) con- 
tains all the interaction terms, and we do not impose any condition like {ipa{k)\ipb{k)) 
or 1. Obviously we have Ea{k), ha{k), Vo(k), hi,{k), V{,(k) and Eb{k) as eigenval- 
ues. Degeneracy due to energy level crossing in k-space implies Ea{ki) = Eij{ki) at 
a certain k point (ki), and at other k points, they are not degenerate. If they are 
always degenerate, then ha{k) = /ib(k), VQ(k) = V(,(k) and Ea{k) = £'b(k), which 
physically mean (pa(k) = ipb(k): this strictly implies that there are two electrons 
occupying the degenerate energy level throughout the k-space (Uh = and ^ = 0). 
As a consequence, degeneracy due to energy level crossing requires /la(ki) 7^ /ib(ki), 
Va(ki) 7^ Vb(ki) and £'a(ki) = Ei,{ki), i.e., Uh = and C 7^ 0. The above logical 
exposition originates from Ref. [55]. 

2.1. Green function formalism. The main feature of the Green function for- 
malism is due to its suitability to describe electronic properties of the interacting 
electrons [IHl HSl 130] ■ The Green function, 

(2.2) g(^t-t') = -'-eit-t'){{c^it),ciit')}), 

defines the probability amplitude of finding a particle at a later time, t in the 
state, a = nkcr, in which the particle was actually created at an initial time, i' 
where t > t' , il is the principal quantum number and a denotes spin-up {(J,a' ,'f,+-^) 
or down (— cr,— cr',J,,— ^). Here, Ca{t) and cjj(t) denote the usual annihilation and 
creation operators for fermions, while the curly bracket denotes anti-commutator. 
Moreover, 

(2.3) e{t-t') = 1 5{t-t')dt. 



Hence, the Green function given in Eq. ( 2.2 ) captures the propagation of an electron 



during the time interval from t' to t [2Hj. By taking the time derivative of Eq. (2.2 1, 
one can show that i2Hl 



(2.4) *^^^^^ = ^ii n 'fit t'){{[co.it),H],cl{t')}), 
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where H is the many-body Haniihonian. If H is the free-electron (non-interacting) 
Hamihonian, then the Fourier transformed Green function [also known as the free 
propagator with superscript (0)] is given by [2511^ 



where, 5 is a convergence parameter. Equation ( 2.5 ) has a pole when E—Ea+iS = 
and the density of states is ^^ S{E — Ea), where S{E — £'„) is a delta function. If 



the interaction is included, say, the e-e interaction, then Eq. (2.5) can be written 

as [2H1ES1 



(2.6) 



GiE,j:) 



1 



E-lEr, 



i6 



T,{a,E)] 

Therefore, there is an additional energy term in the presence of e-e interaction that 
changes the pole structure (spectra), which is known as the self-energy correction 
[I](a,-E)]. In principle, Eq. ( |2.6[ ) captures the overall structure of the solution with 
many-body interaction j28[ 129) . Importantly, Eq. ( |2.5[ ) (without interaction) and 
Eq. (2.6 ) are identical mathematically with Ea + I^{a, E) as the effective many-body 



electronic energy. 

2.2. Hubbard model. The Hubbard model [HI [17] is unique in a sense that it 
was the first model used to interpret the insulating property in transition metal 
oxides. The Hubbard model Hamiltonian [2S] 

(2.7) if« = i?(0)^4,CR, + t ^ c^+d^CR^+UH^nRtriR;. 



R,r 



R+d.CT 



R 



Here, tt-r^ and tir^ are the number operators with spin-up and spin-down, respec- 
tively, i?'"-' is the non-interacting atomic energy level, R is the lattice-point coordi- 
nate, and d is the distance between two nearest ions with certain spin configuration. 
The hopping or transfer matrix elements between the two nearest neighboring sites 
(R,R-hd) is HH] 



(2.8) 



t = / d^r0*(r)v(r)</)(r - d). 



where ^(r) and (j)(r—d) denote the spin-independent (neglected spin-orbit coupling) 
atomic-like orbitals at sites R and R -f d, respectively, which are distorted in the 

crystalline environment, and v(r) = (e^/47reo)(l/r + l/(|r — d|) is the Coulomb 
potential due to those two lattice sites (See Fig. [TV A)). 




Figure 1. (A) and (B) depict the coordinates used to write 



Eqs. (2.8 1 and (2.9), respectively. 
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The interaction matrix elements |28j 

(2.9) Vnn'^l /0*(ri-R)0*(r2-R)-^^^^^^i^0(r2-R')0(ri-R'), 
2 J 47reo|ri-r2| 

gives the Mott-Hubbard gap, Uh, and the coordinates are shown in Fig. fltB). In 



writing Eq. (2.9), we considered the Hamiltonian that describes the two static ions 
and the two-electron interaction, which is a reasonable approximation when one 
only considers the nearest neighbors [28. In this case, the interaction Hamilton- 



ian [ig (from Eqs. |2J| and (2.9|) 



(2.10) -ffilit ~ X! '^RR'4ia4i'a'^R''^"^R'^= X! "^'rR'"-R<t"-rv'. 

R,,(T,R,',cr' R,,cr,R,',cr' 

We have used riRo- = c^^c^a, where Rcr 7^ R'cr', which means c^^c^'a' ~ 0. If 
cr = t, cr' ^ — (t' = I and R = R', then electron rj (electr on at ri with spin down) 
and rl can occupy the same site R, which then leads Eq. (2.9) to be written as 

(2.11) 5I^RR'=E^i^i^^Uh. 

r\,,R, R. 



The Green function based on the Hamiltonian given in Eq. (2.9), in the atomic 
limit is [25] 

(2.12) ^(i?,UH' '-^''-^ ^"-"^ 



E-E(^)+i5 E-[E(0} +Vii]+iS' 



Equation ( 2.12 ) gives two poles where the first pole at E^°^ is due to E—E'-^'> +i6 = 



and the second pole at E^^^i + Uh is from E - [£:(°) -I- Uh] +i5 = 0. 

2.3. lET and the Mott-Hubbard metal-insulator transition. The advan- 
tages of invoking the Hubbard model within the ionization energy theory are four- 
fold, namely, 

(i) We will understand that the screened Coulomb potential for a two-electron sys- 
tem, Vsc (Eq. (19) of Ref. [3T]) captures the same effect as the Mott-Hubbard gap, 
Uh , and it is in the form of Shankar-type screened Coulomb potential [32) , which is 

(2.13) V;e = , ^ ^ e-M(n+..)e^<-" 

where /ies^'"'?) denotes the screening parameter, /x is the constant of proportional- 
ity and |ri — r2| = \Jr'l + r| — 2rir2 cos(6'2). Here, A — 12aB7reo/e^, as is the Bohr 
radius of atomic hydrogen, e and cq denote the electron charge and the permittivity 
of free space, respectively. 

(ii) Since the (V^c) is also a function of the ionization energy (^), one has the ad- 
vantage of studying the evolution of (Ksc) with respect to doping: strength of the 
metal-insulator transition due to doping or changing interaction strength, which is 
not straightforward with Uh. 

(iii) To obtain the Green function as a function of the ionization energy, and 
(iv) New conditions can be invoked such that Uh = and ^ 7^ that do not satisfy 
the free-electron metals. 

For an apparent reason, we have identified ^ interchangeably with the ioniza- 
tion energy and the energy-level spacing. The reason is that the existence of this 
energy-level spacing originated from the screening of the valence electron due to 



6 ANDREW DAS ARULSAMY 

some Coulomb repulsion from an inner electron, which can be observed in any 
two-electron atom or ions, namely, in He, Li+, Be^"*", and so on. For example, the 
Coulomb attraction between a valence electron and its nucleus will increase sys- 
tematically from He to Li+ and to Be^^, which is why we need larger ionization 
energies (due to increasing Z) to remove them. At the same time, we also have a 
stronger Coulomb repulsion between those two electrons as a result of increasing Z, 
which have given rise to a larger energy-level spacing. Increasing Z means increas- 
ing Coulomb attraction between an electron and its nucleus. If Uh = and ^ = 0, 
then the system is in a free-electron metallic phase. The new condition, Uh = 
and ^ ^ gives rise to the strange metallic phase in the normal state of cuprates. 
To see the similarity between Uh and ^ clearly, let us recall the interaction part 



Eq. (2.91 using Eq. (2.131 to obtain 



of the Hubbard Hamiltonian [Eqs. (2.9 1 and (2.10l] here, in which, we can rewrite 



(2.14) {V^n') = ^ y d3rid3r20*(ri - R)0*(r2 - R)[V;c]0(r2 - R')0(ri - R'), 



where the atomic orbitals (0 ) are exactly the same for both Vrr' in Eq. (2.9) 



(undoped) and V^^, in Eq. (2.141 (doped). Using Eq. (2.141, the renormalized 
Green function within lET is 

Therefore, we can now study the evolution of the Mott-Hubbard gap (the shifting 
of the poles) for various doping elements because (1r_r) is a function of ^, which in 
turn uniquely related to the type of elements in a given solid (non- free-electrons). 
In summary, the Hubbard model that gives rise to Uh as a result of e-e Coulomb 
repulsion has been recaptured within the lET by identifying that Uh is equivalent 
to (14c) in the atomic limit. 

Figure [2] schematically captures the similarity between the Hubbard model and 
the lET. For example, for two-ion (A and B) and two valence electron system, the 
Hubbard model needs the accurate atomic orbitals for each electron to determine 
Uh. On the other hand, the lET captures the same scenario via (Ip^), and on 
top of that, the evolution of (Vp^) is also predictable for different ions substituting 
A and/or B. In other words, if we have a system, ABi-^jC^:, then one takes the 
average of ^ as follows; ^ra = {^~x)^b, + {x)£,c , which will evolve with the content of 
ion C that systematically substitutes B. As for the transfer matrix elements, t [see 



Eq. ( 2.8 )] , one does not need to evaluate this parameter since it is straightforward to 
notice the relation, t ex 1/Uh oc l/(T^p) because larger ^r2 implies more localized 
r2 electron and the overlapping between electron r2 with ri is reduced, as it should 
be. 

To evaluate t one needs to know the accurate atomic orbitals [^(ri,R) and 
(?!)(r2,R)] and fortunately, lET does not require such information to predict the 
evolution of (l^pp) for various doping elements. Apart from that, the two solid lines 
(not the arrows or dashed lines) drawn in the lET diagram (right-hand side) in Fig.^ 
is the renormalized interaction of four solid lines in the Hubbard model diagram 
(left-hand side) through (V^p). For example, for the two-ion and two-electron 
system in Fig. [2J electron ee screens electron ca, which is attracted (Coulomb) 
toward ion B, and the strength of this screening depends on (V^pp). Meaning, 
electron Sa can be thought of as an electron brought from say, infinity toward ion B, 
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Hubbard Model 

1/ 





Figure 2. Schematic diagrams for the two-particle Hubbard 
model (left-hand side) and the lET version of the Hubbard model 
(right-hand side). A and B are two different static ions with one 
valence electron each. The separation between the electrons is 
|ri — r2|, and R is the separation between the ions. The effective 
charge, e of the electrons in the Hubbard model depends on the 
accurate atomic orbitals [^(ri,R) and ^(r2,R)], while in the lET 
version, the effective charge, e is captured by the averaged screened 
Coulomb potential through the parameter, ^. 

and its attraction is determined by (Kc'^)- Here, the two solid lines are captured by 
(Kc'^)- Simply put, electron ca is attracted toward ion B and screened by electron 
cb, while at the same time, electron cb is attracted toward ion A and screened by 
electron ca, which then determine the Coulomb repulsion between electrons ca and 



cb- As a consequence, using Eq. (2.13) (from Ref. [SI]), we obtain 



o2 



(2.16) Kt = 1 r -e-A"('-i+'-^)^"'''^"' 

47reo|ri-r2| 

: Bb is attracted toward A and screened by ba. 



-4a({b) 



(2.17) V:^ = ^ -^.{r^+r,)e 

47reo|ri-r2| 
: ba is attracted toward B and screened by ee. 

Here, ^a and ^b capture the attractions between electrons A and B toward ions A 
and B, respectively. If ^b > £,A then V^ > V^^, which implies that electron A is 
more screened (less localized) than electron B (this statement has been proven in 
Ref. [3T] for atomic He) that eventually leads us to Fig. [SJ 

FigurelsFa) tells us the degree of overlapping (shaded section) between the atomic 
orbitals A and B. The atomic orbital, A is less localized and therefore, is widely 
spread as compared to the atomic orbital, B. Figure [3Fb) is for the hypothetical 
system ABi-^^C^^. If ^a < £,c < £.b, then in average, the shaded area (degree of 
overlapping) increases with increasing content of ion C. Add to that, the mean 
positions of the shaded area should be nearer to the ion B (due to strong Coulomb 
attraction between electron A and ion B). Note here that, ^b > Ca also implies 
strong Coulomb repulsion between the attracted electron A toward ion B and elec- 
tron B. Therefore, in the lET framework, the Hubbard Coulomb repulsion is given 

by 

(2.18) UH-(yAR)«(v;c> + (Kc>- 

Importantly, if for example, ^a — ^ oo, ^b — >■ oo and ^c ~^ oo, then the respective 
ion cores (A, B and C) are rigid with an effective positive charge. If on the other 
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(r) (a) 



<^(r)t (b) 




Figure 3. Schematic diagrams for the overlapping of atomic or- 
bitals, A and B. (a) The shaded area, which represents the degree 
of overlapping depends on the relative magnitudes of ^a and ^b, 
where ^b > £,a- (b) If the conditions, ^b > £,A and ^c < £.b are sat- 
isfied, then the degree of overlapping for the orbitals A and (B,C) 
increases with the content of ion C. Here, B,C implies Jii^^Cx- 



hand, ^a ^ 0, £,b ~^ and ^c — ^ 0, then Eq. (2.18) satisfies the Thomas-Fermi 
approximation (screening of free-electrons) . 



Unfortunately, Eq. (2.15), even though captures the evolution of the poles with 



changing interaction strength or energy level spacing, it is still physically identical 
to Eq. (2.12) in which, Uh = also implies (V^'^j^) = regardless whether ^ = or 
^ 7^ 0. Therefore, we need to find an alternative way to incorporate Uh = and 
^ 7^ into the Green function. A new direction is possible if one employs the total 
energy eigenvalue given in the lET-Schrodinger equation [3T], Hip — {Eq ± ^)lp, to 
simplify the Hubbard Hamiltonian (Eq. (|2.7[)) such that 



(2.19) 



-"lET 



^'°^ Yl 4iaCR<T + C XI "Rt"R;- 



R.i 



R 



Here, iJ^^-* = Eq where Eq is the total energy for OK. The plus sign (in ±^) is 
for the electrons, while the minus belongs to the holes. These signs imply that the 
numbers of electrons and holes are calculated in the limits to (X) (/„ ) and — cx) to 

(/_ ) , respectively. Further, we have conveniently and correctly ignored the term. 



tE 



R-l-d,CT ^R+d(T 



cro- because t is redundant within the lET formalism in which. 



^ is the responsible parameter to take into account any change in the excitation 
or hopping probability. Having said that, we can derive the more user- friendly 
renormalized Green function 



(2.20) 



GiE,0 = 



1- n_ 



E~EW+iS E-[E(o') +^]+i6' 



This time we have to look up the Appendix for a detailed derivation (of Eq. ( 2.20 )) 



because of irritating algebra. In any case. Fig. [4] indicate the predicted density of 
states when ^ > 0, in which, the new condition Uh = and ^ 7^ is actually related 



to Eq. (2.1), and is expected to be responsible for the strange T-linear metallic 



property of cuprates above Tgc. The full implication of Eq. (2.20) in this respect 
are exposed in the following Section. 
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Figure 4. Predicted density of states, D{a,E) based on the 



atomic- limit Hubbard model within lET or Eq. (2.20). Diagram 
(A) depicts the usual insulating behavior similar to Uh where ^ ^ 
and Uh 7^ 0. In this case, both ^ and Uh represent the band gap. 
Whereas, diagram (B) shows the Fermi gas or free-electron metal- 
lic behavior because ^ = Uh = 0, in other words, the band overlaps 
due to energy-level crossings. However, strange metallic behavior 
is anticipated when the energy levels cross (or the energy bands 
overlap) in such a way that Uh — 0, while ^ remains finite (^ 7^ 0) 
as indicated in diagrams (C) and (D). 



By now, we may have noticed that within lET, the atomic ionization energy, or 
the atomic energy-level spacing replaces the role of wave function (required in the 
Schrodinger representation), or the role of Green function in the Green function 
formalism. However, there is an intrinsic inadequacy within lET because it cannot 
be used to calculate the energy-level spacing. Instead, lET exploits the atomic en- 
ergy levels to predict numerous physico-chemical-biological processes for molecules 
and solids via the energy-level spacing renormalization method. This means that, 
one still needs to make use of the Schrodinger equation or the Heisenberg repre- 
sentation, or the Green function formalism to calculate the energy-level spacing of 
any quantum system. 

3. Normal-state properties of cuprate superconductors 

There are numerous proposals that are available to explain the origin of T- 
linear property and superconductivity in solids, namely, the resonating valence 
bonds [31 [Ml [33 ES], the pure and marginal Fermi liquid theory [5711551 15^ 1 liO] . 
the non-Fermi liquid phase UH [H [IS [501 SI [13] , the hidden Fermi liquid with 
different scattering rates HH [13 US] , and the quantum critical-point phenomenon 
and the Hubbard model [HI [HI [47l [48l [H [5^ , Quite recently, we also have the 
topological superconductors (SH [52l [53]. In this work however, we focus on lET 
and find out how to exploit this theory to explain the well-established experimen- 
tal results self-consistently, and accurately. Recall that we will ignore the origin 
of (i) superconductivity, (ii) the pseudogap phenomenon, and (iii) the change of 
resistivity with temperature, giving rise to the T-linear effect for high T, and the 
approximately T^-effect for very low T under high magnetic fields. 

We deliberately ignore the origin of T-linear effect in the normal state of super- 
conductors, not because we are unable to provide any explanation to its existence 
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(in fact, we did provide in Refs. [Ml [55l |56] ) , but because the origin of T-linear effect 
cannot be pinned down to a single cause. The obvious cause can be due to different 
scattering rates (l/Tt,ansport oc T and l/Tfjaii oc T^) as proposed in the hidden Fermi 
hquid theory or due to boson-boson interaction with l/rboson oc T proposed by the 
gang of Phihps (Phihp Anderson-Philip Casey |45l|46] and Philip Phillips [47ll48] . 
respectively). Apart from the Philips, it was Kastrinakis who had proposed a T- 
linear effect in cuprates (both above and below Tgc) entirely from the interacting 
Fermi liquid theory in the presence of the van Hove singularities [371 I3H]- The 
other reason is due to the change in the valence states of multi-valent element with 
decreasing temperature. Changing valence states means changing ^. Experimen- 
tal proofs for this second reason were reported by Dionicio [57] using the results 
of Fukuda et al. [SB], and the data points have been reproduced here in Fig. [5] 
Such a dependence gives rise to a T-dependent carrier density (see Eq. (3.1|), and 



together with different scattering rates, they influence the resistivity versus tem- 
perature curves. For example, if the carrier density, n ex T, while l/Tdoctron oc T^, 
then we have the famous T-linear effect [SH [55], p{T) ex T. Unfortunately, one 
also need to explain the Hall angle that is proportional to T'^, as well as the almost 
temperature-independent heat conductivity [SHI 1101 IM] self-consistently. 
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Figure 5. Temperature-dependent valence state of the element 
Eu in EuCu2(Gei_a;Si2;)2 due to increasing x, which were inferred 
from the Eu-Lm edge X-ray absorption spectra obtained from 
Fukuda et al. [58] . This figure was carbon-copied from the work of 
Dionicio [57] . 



What we attempt to do here are to explain the existence of thermal assisted 
^QPT in cuprates above Tg^ and to find out how does this continuous QPT is 
responsible for the existence of non-metallic (insulating), strange metallic, Fermi 
liquid (go > 0) and free-electronic (Fermi gas; go = 0) phases. Here, go denotes the 
bare coupling constant following Ref. [3], and we will get to this point in detail later. 
First, we need an appropriate figure to understand the existence of these phases, 
and Fig. l6lA) depicts the relevant phase diagram for cuprates drawn schematically 
following Refs. [Ml ISH IMl IMHSSISS IS7] . Figure [6];B) shows the resistivity versus 
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temperature {p{T)) curves for the paths, II to V drawn in (A). The dashed hues in 
Fig. IgI^B) show the deviations of p{T) from the T-linear curves for paths II and III 
at Tj*! and Tj*jj, respectively. As stated earher, we will ignore any p(T) curve that 
contain the pseudogap phase between Tgc and T* because such a phase is never a 
prerequisite for superconductivity in cuprates. 
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Figure 6. (A) Schematic phase diagram for cuprate superconduc- 
tors plotted as a function of the ionization energy where the free- 
electron system is uniquely defined when ^ = 0. If one increases 
^, the system may acquire interacting Fermi liquid phase, followed 
by the strange-metal and pseudogap phases, before achieving the 
insulating antiferromagnets, through a series of electronic phase 
transitions at a constant temperature. (B) Resistivity (p) versus 
T curves follow the paths, II, III, IV and V as sketched in (A) 
in which, these schematic curves have been observed experimen- 
tally [SHI IMl ISH ISSl ESI EZ] . All the arrows pointing downward in 
(B) indicate Tgc, while the upward arrows denote T* . The metal- 
to-insulator crossover temperature is labeled as T„- 



To keep the discussion straightforward, we need to exploit pi^) curves along the 
paths, IV and V that may also include the metal-to-insulator crossover temperature 
(Tcr). In addition, such curves should also be presented such that TgcS have been 
suppressed by means of magnetic fields so that the dotted lines in Fig. IgIB) can be 
confirmed. These dotted lines do not shoot to infinity when T — >■ 0, instead, they 
approach finite values that correspond to a metallic (Fermi-gas- or Fermi-liquid-like) 
behavior. The experimental data points satisfying paths IV and V in Fig. IslA) are 
depicted in Fig. ^ which were measured by Takagi et al. [66 , Cooper et al. 67] and 
Hussey [5H1- These p{T, x) curves are for La2-xSra;Cu04 polycrystals (A) and single 
crystals (B). There are two distinct electronic properties that need to be tackled 
self-consistently, (i) the downward shift in p{T, x) due to increasing x and (ii) the 
transition to the Fermi-liquid-likc metallic behavior after the T-linear effect in the 
absence of Td-. The answer to (i) is unambiguous and self-consistent [5^ IB51I7D] . 
For example, the change in the carrier density 



(3.1) 



(r,0= / f{Eo,OD{Eo)dE = 



knT 






exp 



Scuprate 

keT 



is responsible for the downward shift of p{T) with increasing x (see Fig. u\ because 
^Sr2+(807 kJmol"^) < ,^La3+(1152 kJmol^^). Here, ^cuprate is now the energy level 
spacing in La2_a;Sr2,Cu04 determined by the content of Sr and La, hence Scuprate = 
a;^(Sr^+) + (2 — a;)^(La^"'") where Scuprate is proportional to their atomic ionization 
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energies of Sr^"*" and La'^^. Moreover, /(-Eo, i^ the ionization energy based Fermi- 
Dirac statistics given in Ref. 31 , D(Eq) = ml/Trh? is the two-dimensional density 
of states, ml is the electron effective mass, and Ep is just a constant and we take 
it to be zero for convenience. 

We show here that there exist a ^^QPT from the strange- to Fermi-metal (with 
finite conductivity when T — >■ 0) within JET. Recall here that lET cannot be applied 
to Fermi liquid or Fermi gas because the energy-level spacing, ^ is either zero, or 
does not play any significant role for such systems. When one looks closely, the 
p{T, x) curves in Fig. (TT A) and (B) cross each other, violating the prediction from 
Eq. (3.11 for different x and for T < lOOK. The first reason that comes to mind 



is the T-dependent scattering rates becoming the dominant mechanism as a result 
of ^ now has been reduced to an insignificant constant. Note here that any Fermi 
liquid or Fermi gas system has the power to render ^ useless and/or redundant 
as it is suspected to be the case here for T < lOOK. Nicely enough, this violation 
exists without significant changes occurring to ^. For instance, when T > lOOK, 
^ very accurately determines the shift of the resistivity curves with x, however, 
for T < lOOK, the effect of ^ on p{T, x) is strangely switched-off (not because 
^ — >■ as r — T' OK). For example, the above transition for p(T) from strange- 
to Fermi-metal does not require significant changes to ^, instead ^ continuously 
ceases to be relevant. This means that the Fermi metal and its scattering rates 
become increasingly significant for T < lOOK without requiring ^ — > 0, and this 
conclusion somewhat points in the direction of the hidden Fermi liquid theory 
proposed by Anderson and Casey [451 146] , as well as the interacting Fermi liquid 
theory advocated by Kastrinakis [371 138) . Unfortunately, we are unable to pin point 
the transition temperature between the strange- and the Fermi-metal, Tp^ where 
Tp^ is definitely not related to lOOK discussed above. In fact, the value lOOK is 
the temperature when p(T, x ~ 0.23) curve crosses the curve for p{T, x = 0.26) (see 
Fig-^B)). 
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Figure 7. (A) The resistivity versus temperature (p(T, x)) curves 
for the polycrystalline La2-a;Sra;Cu04. (B) pa\^{T,x) curves in the 
ab-plane (Cu02 plane) for the La2_a;Sr2;Cu04 single crystals mea- 
sured at high magnetic fields to suppress Tgc- Both (A) and (B) 
figures were carbon-copied from Ref. |45j in which, the experimen- 
tal data points were obtained from Refs. [551167) . 
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Anyway, this so-called ^QPT (from strange- to Fermi-metal) is observable in 
the p(T,x) measurements (see Fig. ItFB)) because p{T) — ^ finite value as T ^- 0. 
Without this transition, p{T) — >■ oo as T — )■ due to the existence of ^ 7^ in 
Eq. (3.1 1. This §^QPT is actually an electronic thermal phase transition due to the 



existence of Eq. (2.1). For example, for a two level system, a single energy level 
crossing at ki technically implies ^(ki) — |ft,a(ki) — /ib(ki)| -f |va(ki) — v;,(ki)| 7^0 
where Eai}f.i) = £'fc(ki) and Uh = by definition. Applying this knowledge to the 
above problem means, for T > TpJ^, the electronic conductivity requires transport 
from one quantized energy level to another, which needs to overcome ^(k) at the 
points where the energy levels cross. This gives rise to the strange metallic behavior 
in accordance with Fig. WTC and D). Whereas, for T < TpMi the same electronic 
transport from a quantized energy level does not cross to its neighboring energy level 
at the point of crossing. As such, ^(k) does not play a significant role below T^^, 
giving rise to the Fermi-liquid (not Fermi gas) metal, and its scattering rates are the 
ones that determine the temperature-dependence of p{T). This low temperature 
Fermi metal is a Fermi-liquid type, rather than the Fermi-gas (free-electron) type 
because the latter strictly and literally requires ^ = 0, even though both Fermi 
gas and Fermi liquid have continuous network of Fermi surface. A continuous 
Fermi surface exists only for 1^ = and/or when 1^ 7^ and T < Tp^ where TpJ^ 
denotes the ^^QPT temperature between the strange- and the Fermi-metal. Any 
other condition does not allow a network of continuous Fermi surface independent 
of (,. Apart from that, this thermal- assisted quantum phase transition includes the 
continuous (higher-order) thermal phase transition, and consequently, Tp^ cannot 
occur at a fixed temperature, instead, it represents a temperature range without 
any discontinuity that can be observed indirectly in p(T) measurements (when 
p{T) — > finite value as T — > OK). Therefore, we can claim here that all second- and 
higher-order TPT are thermal-assisted ^QPT. On the other hand, the claim that 
all first-order TPT must go through a thermal-assisted '^QPT at a constant T has 
been proven elsewhere [7T| . 

There is however, one other issue that escaped our analysis thus far— the origin of 
Tcr depicted in Fig. ptB) for path II indicates an insulating behavior just above the 
superconducting transition temperature. This non-metallic behavior is observable 
because ^ > /csTcr where ^ 7^ here does not imply the existence of a band gap. In 
the presence of a well separated conduction and valence bands, and in the absence 
of any energy level crossing between them, ^ 7^ is the energy level spacing related 
to a band or energy gap (Eg). Therefore, for any ^ such that ^ 7^ and Eg 7^ 0, one 
will always obtain an insulating system, whereas, for ^ 7^ and Eg = 0, one instead 
captures the strange metallic behavior at high temperatures, and possibly a Fermi- 
liquid metal at low temperatures if and only if there is a ^QPT. The existence of 
Fermi liquid and strange metal for different temperature ranges implies that the 
emergent phenomenon put forth by Anderson [27] is actually related to QPT. In 
other words, indeed "more is different" but we can certainly add details to it, and 
claim, more interactions (not necessarily competing) within a given system will lead 
to different types of quantum phase transitions and emergent phases. Finally, path 
I in Fig. pTB) shows the quantum phase transition due to doping or changing ^, 
and in this case, the '-'QPT is from an antiferromagnetic insulator (large ^) to a 
strange metal (low ^) for T > T* . 
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Having said that, we continue and complete our analysis on the normal state 
properties of cuprates with respect to the proof developed by the trio, Parameswaran, 
Shankar and Sondhi [3^ . We recall the definition of Fermi-liquid based on the bare 
coupling constant (go), and the s-wave Cooper-pair frequency (fi) derived by them. 
The relevant result is the two-point Cooper-pair correlation function [3], 



(3.2) r(fi;go,Ashankar) = ;^ 

ZTTVp 



a\ log (f^/Ashankar) 



l+ago\ log (fi/Ashankar) I . 

where fl = (fi^ -I- P^)^/^, P is the s-wave Cooper-pair momentum, a denotes a pos- 
itive constant, Ashankar IS the Shankar cutoff parameter, Vp denotes the Fermi 



velocity, and Eq. (3.2 1 must satisfy Q < Ashankar- For Fermi gas (free elec- 
trons), go = and for T = OK, one has r(f7; go, Ashankar) « | log (i7/Ashankar)|, 
which diverges logarithmically as 51 — >■ 0. On the contrary, for Fermi liquid, 
r(fi; go, Ashankar) -^ l/27rVF5o, which converges to a constant as fj — > 0. This result 
clearly indicate that Fermi gas is indeed a special case within the Fermi liquid for- 
malism such that when ^ has been switched-off, one can simply take ^ = because 
^ is now an irrelevant constant. Obviously, there is a one-to-one correspondence 
between the Fermi gas and Fermi liquid when ^ — >■ 0, which is also in agreement 
with our earlier claim [S2] that the Fermi gas (or a free-electron metal) has got to 
be an emergent phase. For example, we can write [5?| Ashankar = {hk"^/ 2m* )Aiet 
and replacing Ashankar with (/ifc^/2m*)AiET in Eq. (3.2) does not change the above 
divergence or convergence in any way. In particular, for the Fermi liquid systems, 
AiET — exp (A^) is a constant because ^ is itself a constant, while A = (127reo/e^)aB 
is just a collection of constants defined earlier. As for the Fermi gas, Ajet — ^ 1 as 
^ — ^ 0. Therefore, replacing Ashankar with {hk'^/ 2m* )Ai-et within a Fermi gas or a 
Fermi liquid system does not change anything physically, as it should be [SU |S2] , 
except that Fermi gas is a special case exists when and only when ^ = 0. 

4. Further analysis 

Contrary to our models derived earlier [21] [22l |24j , the particles responsible for 
the strange metallic behavior discussed in this work are still electrons, but with 
strong correlation such that it requires ^ 7^ and ^ is not redundant. Free-electron 
metals simply require ^ = 0, while Fermi liquid embeds the condition, C 7^ but 
its magnitude is an irrelevant constant. Therefore, our strange metal is due to 
strongly correlated electrons requiring ^ 7^ 0, and ^ plays an active role here such 
that the electron-hopping from one energy-level to another needs to overcome ^. 
These energy levels are crossed for strange metals, but if they are gapped, then one 
obtains a semiconducting (or an insulating) normal state, however, this depends 
on the relationship between ^ and T. For example, one gets a metallic property if 
^ < T, and a semiconductor otherwise. 



The Green function in the atomic limit given in Eq. (2.201 is not applicable for 
Mott insulators. This is a false assumption. Proving it false is the main aim of this 
final piece. We will proceed to demonstrate why and how ^ captures the hopping of 
particles with different elements and elemental composition, without violating the 
particle-hole symmetry, without any ad-hoc physical interpretations. Our analyses 



follow Ref. [72]. We show Eq. (2.20) correctly exposes the spectral weight transfer 



mechanism in accordance with the data obtained from the soft X-ray absorption 



spectra, measured by Chen et al. (73] ■ In fact, in the atomic limit, Eq. ( 2.20 1 requires 
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Figure 8. Soft X-ray absorption spectra at the oxygen-K edge 
for the sample, Lajt^Sr^+Cu^t^Cu^+O^^^j obtained by Chen et 
al. |73| . The main figure with experimental data points was carbon- 
copied from Ref. [73] . The intensities (/) for these peaks have been 
normalized, backgrounds subtracted, and the data were fitted with 
a Gaussian function. There are two main peaks, one being the 
high-energy (HE) peak (occurring between 530 and 531 eV), while 
the one to the left is the low-energy (LE) peak occurring slightly 
below 529 eV. The photon energies labeled with p and q show the 
red-shift (from p to q) in the LE peak for increasing x, while the 
blue-shift occurs from r to s in the HE peak for the same doping 
concentration (0 < x < 0.15). These shifts satisfy the inequality, 



p > q > r > s where q — p — A 



red 

shift ' 



= A 



blu 
shift ' 



9 = A. 



and s — p = A„iax. Here, A^J^ff^ and A^^^f^ denote the total red 
and blue shifts for the stated doping levels, while Amjn and A^ax 
are the minimum and maximum changes to the electrons energy 
levels, again within the investigated doping range. The change in 
the normalized intensities for the high- and low-energy sectors are 
denoted by A/he and A/le, respectively. 



one to write ^ = CI'ysTms^ for molecules, one has Q°J^™i^' 



namely, for cuprates one writes ^cuprate 



and for a specific solid. 
Hence, the many-body problem here is 



to calculate the many-body energy level spacing denoted by Cb^^J- For atomic 
systems, ^ is known almost exactly since we can obtain it from the atomic spectra, 
whereas, for other compounds and molecules, we need to invoke the ionization 
energy approximation [54], which states, the energy level spacing of a compound 
is proportional to the energy level spacings of the compound's constituent atoms. 
For example, the energy level spacing for a hypothetical compound, ABi_a;Ca; is 
given by ^aBi_^c^ = x^c + (1 — x)£,b (see Fig. [3J(b)). In this case, we can exclude 
the contribution from A because ^a is a constant anyway. 

We now recall the soft X-ray absorption spectra measured by Chen et al. [73] 
depicted in Fig. [HJ These spectra correspond to the changes in the concentrations 
of Sr2+, La3+, Cu2+, Cu3+ and O^" in the compound, La^l^Sr^+CuJlj^Cu^+O^;^. 
Unlike the simplified hypothetical ABi.j-C^, system discussed earlier (the concen- 
tration of A is made to be constant), doping Sr^+ into Laj^Cu^l CUy+04"^Q5 gives 
rise to changing 6 and y, which have been discussed in detail elsewhere for other 
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oxides and mult i- element solids [70l [74l [751 [76] . For example, such changes are due 
to defects, types of defects and the required overall charge neutrality in the system. 
In particular, since 5 ~ 0.005, we need y = 0.01/3 for charge neutrality. In the 
presence of Sr (0 < x < 0.15), both y and 6 can, and will change, satisfying the 
linear equation, x — y ~ 8 — 2S where {x, y} e [0, 1] and 5 £ [0, |]. This knowledge 
will become handy shortly. 

Now, one can immediately observe four interesting features in the spectra plotted 
in Fig. [Hj The first being p and q for x > appear at energies lower than r and 
s where p — q — Oiix = (undoped sample). Second, the normalized intensities, 
/le increases, while Jre decreases with x. Third, the red (Agf^'j'fj) and blue (A^j'^^j'^.) 
shifts appear for the low- and high-energy peaks, respectively, with increasing x, 
and finally there are some asymmetric changes to the normalized intensities for the 
range < a; < 0.15, such that A/le > A/he (see Fig. [8|. In other words, /le 
increases faster than jd/HE/da^l- 

The answer to the first feature is straightforward, it is due to the fact that ^31-2+ 
(807 kJmol"^) < £,i^a_3+ (1152 kJmol"^). This means that the electrons with low 
ionization energies (from Sr^+) need low photon energies to be excited, compared 
to the electrons with large ionization energies (from La'^"''). Here, all ionization 
energy values were averaged following Ref. |54j . while the values before averaging 
were obtained from Ref. [77. The second feature is also obvious, the change in the 
normalized intensities are due to increasing x and decreasing 2 — x for the low- and 
high-energy sectors, respectively. 

The red and blue shifts in the spectral weight transfer are also depicted in Fig. [8] 
with appropriate arrows. Understandably, the red-shift originates from the con- 
centration of Sr^+, or due to increasing x. This implies that more electrons from 
Sr^+ are needed to establish this LE peak at the lowest possible energy as a re- 
sult of ^Sr2+ < ^La3+ • Rccall that lET requires electrons with the lowest atomic 
ionization energies to form the valence electrons in a non-Fermi gas compound con- 
taining that particular atom. These (valence) electrons also interact weakly with 
the core electrons (coming from La and Cu), giving rise to smaller energy level 
spacings. The core electrons form large energy level spacings due to strong interac- 
tion among the core electrons [31]. As a consequence, the HE peak corresponds to 
the core electrons (due to large ionization energy) coming from La2_^. These core 
electrons interact more strongly with other core electrons coming from CU]^_ and 
Cu^+ (Ccu2+ = 1352 kJmol-i and Ccu3+ = 2086 kJmor^) giving rise to HE peak. 
This strong electron-electron interaction also explains why the energy distribution 
can and will spread from low to high energies as a result of the interaction among 
the core electrons and between the valence and core electrons. 

We know from the above discussion that increasing x causes y to increase to 
maintain charge neutrality, and therefore, the decreasing number of large ioniza- 
tion energy electrons from La|t^ are systematically being compensated by the large 
ionization energy electrons from Cu'^"'" (due to increasing y) . This scenario is further 
enhanced if 6 is also found to increase with increasing x. Therefore, we can now see 
the reason why HE peak blue shifts and A/le > A/he; they are due to increasing 
y, 2 — X > 1 and/or increasing S. For example, in the absence of this compensation 
effect, the HE peak should red shift with decreasing 2 — x. But this is not the 
case, made sure by increasing y, and the peak will not disappear even when x = 1 
because 2 — a; 7^ and y oc x. Note here that we did not assume La^+ contribute 
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three electrons, while each Sr^+ gives two electrons. As a matter of fact, the carrier 



density increases with increasing Sr^+ in accordance with Eq. (3.1). Therefore, the 
compensation effect slows down /he from decreasing faster than the increasing /le, 
and therefore A/le > A/he- Indeed, the blue-shift is due to ^cu2+.3+ (1352, 2086 
kJmol^^) > ^La3+ (1152 kJmol^^). In summary, we did not encounter any tragedy 
between lET Green function and the spectral weight transfer measurements, in- 
cluding the observed doping-dependent spectral red and blue shifts in the so-called 
La-Sr-Cu-0 Mott insulator. 

5. Conclusions 

Our primary objective is to prove the existence of a thermal-assisted quantum 
phase transition between a strange- and a Fermi-metal in cuprates. To achieve this, 
we have first developed the proof where an electronic phase transition occurring 
in the normal state of cuprates is a type of ^QPT. This transition is exposed 
to exist in cuprates above the superconducting transition temperature, when the 
temperature is lowered in the presence of high magnetic fields, using the Green 
function formalism. Using the renormalized Green function, and by requiring (i) 
zero energy gap, (ii) non-zero energy level crossings, and (iii) non-zero energy-level 
spacings, we obtained the strange metallic phase with its resistivity curve shifts 
accordingly with respect to changing energy-level spacing or doping. At a lower 
temperature however, there is this ^QPT, which is responsible for the strange- to 
the Fermi-metal transition as a result of insignificant effect from the energy-level 
spacing because it is being switched-off continuously with lowering temperature at 
the points where the energy levels cross. 

Along the way, the proof developed by Parameswaran, Shankar and Sondhi is 
coupled to the ionization energy theory to show that indeed Fermi gas is a special 
case within the energy-level spacing renormalization group method because Fermi 
gas can only exist if the energy-level spacing is zero, while the original condition 
that requires infinite energy-level crossings is actually necessary, but not sufficient 
for this phase (Fermi gas) to exist. 

Finally, we demonstrated why and how the spectral weight transfer occurs for 
different doping elements and concentrations, including the red and blue shifting 
effects, and the asymmetric changes to the normalized intensities between the low- 
and high-energy peaks. We found all the features observed in the soft X-ray spectra 
for the La-Sr-Cu-0 system have been captured self-consistently. These spectra, and 
their red and blue shifts nicely correspond to the changes in the ionization energies 
for different doping elements and concentrations. 

It would be interesting if we could find ways to exploit this energy-level spacing 
such that it can be associated to the origin of superconductivity and the pseudogap 
phase, which have been excluded systematically for convenience, but this is another 
story for another day. 
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6. Appendix: derivation of Eq. (2.20) 



We start from Eq. (2.2 1 



(6.1) 
Subsequently, 



dt 



g{t-t') = --9{t-t'){{c^{t),cUt')}) 



-^^9^(c.(t)cl(t )) - -e{t - 1 ){-^^ci{t ) 



(6.2) 4^^(4(0c«W)4^(i-t')(ct(0%f^). 

Noting Eq. (|2.3|) and after some algebraic rearrangements, 



(6.3) ^/^^^^^ = 5{t t')({c«(t), 4(i')}) - le{t - f){\^^n^-^,ci{t'] 

Now, the expectation value of operator Cq at time ti is 

(6.4) (c„(i)) = {^{t)\c^{h)\^{t)). 
The equation of motion with respect to a constant Ca{ti) is 

(6.5) ^n'^ = {^n'^)c.it,mt)) + (^(.)|c„(.0 (^^^^ 

where ip{t) is some time-dependent wave function and therefore, \^{t)) is the state 
vector. The complex conjugation requires, {ihd{ip{t)\)/dt = {ip{t)\{—H), thus 

.^5(c„(t)) 



dt 



= (vp(i)|(-i?)c„(ii)|^(t)) + {ip{t)\c^ih)H\ip{t)) 



(6.6) 



= {^mCa{tl),HMt)). 



We now replace the dummy time label ii with t and substitute Eq. (6.6) into 

_dg{t-t') 



Eq. (6.3) to obtain 
(6.7) ih- 



dt 



1 



e{t-t')(\[c^{t),H],cUt') 



where H = Hfj^j, from Eq. (2.19), and {cj^(i),c]^(t')} = 1 because 6{t - t') = 1 
We move on by first evaluating 

[c„(i),i:fiET] - c,4^(°)4,c,„- +e4'C,^,c]_,,c,_,,) 



+i,[Cia-C-,Cj 



'-■j<j''--jcr '-■^'^J 



C^-./ C^'fj' C ■__/€,, 



-jrr' -JO Cj-o-'Cj-cr' ^ja-'^jcr' ^j-a'^J-cr' ^^i^ 



)^ 



where i and j are dummy indices (dummies for short) replacing R, while the spin 
dummies, a' and —a' denote spin-up and down, respectively, however, the dummy, 
a can either be spin-up or down. All the annihilation and creation operators are 
summed over and these dummies are necessary to keep track of these operators 
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when they operate from left to right or vice versa. The term with E^^'^ on the 
right-hand side (RHS) of Eq. ^^ is 

(^D.yj -C/ \pj(T' I Cj(jfCj(T'Cicr C.jfj'Cj(j' Ci(j jOijO(j' (J — 11/ Ci(j . 

The first term with S, on the RHS of Eq. ( |6.8[ ) is 



Cjg.: CiaCjcr' Cj_^, Cj^ 



(6.10) 






The above and ah the similar evaluations follow the fermionic-operator rules, for 

example, Cia-C-, l C- /C^^j, CifjCja-' — Cja' ^ia 7 ^ja' — ^ja'^3^' ^^^ "^j — a'^ja' — 

Cja'Tij-a'- Anyway, the last term in Eq. (6.8 1 is nothing but the last term on the 
RHS of Eq. ( [OOl ). Therefore, 



^ja' '^j'^' ^j-a' ^'j-<y' ^ii 



(6.11) 






In the last step, we have replaced the dummies, j, —a' and a' with i, —a and a, 
accordingly via the Kronecker deltas and ni„Ci„ = cl^CiaCia — 0. When the above 
creation and annihilation operators operate from left to right or vice versa, two 
Kronecker deltas pop out in response for each operation, one for the lattice sites 
and the other for spins. We can now write (after replacing a with dummies, j and 



ih 



dg{t - 1') 
at 



l-'-0{t-t'){{{E(^'^c.,,+^rM^,c,,),clit')}) 



= l-i?(")l0(i-i')({Q.(i),4(t')}) 



.12) 



-e-0(t-t')(K_,c,,(i),c]„(t')}). 



Using Eq. (6.1) we can rewrite Eq. (6.12) 

a 



(6.13) 

Note here the new substitution 



at 



(6.14) 



e(l)(t-t') 



g{t-t') = i + ^g'^'>it-t') 



{t~t'){{n,^„c,„{t),cUt')}}, 



which needs to be evaluated in the form of 

'^~dt =^^*' t'){{n^^ac.a it), cl (<')}) 



.15) 
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We have dropped ^ from Eq. (6.15) for the time being. Now, since 5{t — t') = 1, one 
can then employ the fermionic rule, {cicr(i), dg.(t')} = {cia,Cj^} — 1 and Eq. (6.6) 
to obtain 



.16) ^fl ^^^'^^g^ *'^ - {n^-a) - ^^(i-i')( \ K_.C,.(t),ijfET], 4(^0 



The next obvious step is to evaluate 



We first evaluate the term with E'^^^ from Eq. (6.17) 

E^°^n, 



CiaCj^iCja' ^^ja'^jc'^i' 



.18) E'-°^n,_ 



\^ Cjfj'^ia )Cj(T' C ■ / Cjct' Cjcr 



'ja' '^ Cj^iCja'Cia Cj^,Cja-'Ci^ 



S-S / = E^°^n- c 



^l — (7^1(7 • 



The term with ^ from Eq. (6.17) is 



'^i — (T^io'^j(j'^ja'^j — (j'^j — (T' ^ja'^j^' ^j — (j'^j — a''^i—(T^i(T 






J ... J 



J .. J 



Ci(7Cj^'Cjcr'Cj_^,Cj^r7' C.-^iCjcr'Cj_„iCj-cr'Cia 



Cj(j'Cj^„iCj-a' + Cjo-zCjo-'Cj-o-' + C^-g./Cj(j'C^-_o-/Cj_CT'CiCT 



'^ja' '^jf' S-o'' ''J ^'^' ^i' 






'^j — a-'^j — a'^ij^cra' ~r C ■ /Cj^r'^^j — cr' ^ij ^crcr' ^a 



-jcr'^j 



C-ia^i — (j^i — (T ~r C^^CicrCia- 



(6.19) =^ni_crCio-, 

again the Kronccker deltas took care of the dummies, (cia) — and 



T-i-a ^ '^'i-a'^'i-a'^'i-a'^'i-'y ^ '"'i-a\^ C,i_^Ci^a- ) ^i-a 



(6.20) 

Therefore 

(6.21) 



„t J 



[n^-aC^a{t),Hf^rJ,] = in,_aC,„{t) + E'^°'> m^aCiait) , 
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and after substituting Eq. (6.21) into Eq. ( |6.16[ ) 



ih- 



0(t-t')^|[«.-.C.,(f),i/fET],4(^')} 



.22) 



(ri.-.> - e;^e(i - 0( I N-^c,,(t), c]^(i') 



-i?(")-0(i-<')(^N-.c,.(i),4(i') 



Using Eq. (6.14) one can simplify Eq. (6.22) 



at 



.23) 



*s-^'°'-« 



5<«(( -(■) = {'•.-,)- 



Subsequently, one obtains (after using the correspondence rule, ihd/dt — > E) 



.24) 



g(^\t-t') = 



E-E(") -f 
We substitute Eq. (|6.24) back into Eq. (|6.13l to arrive at our final destination 



(i?-i?("))g(t-0 = l + ^%^. 



(i?-i?(°))g(i-i') = 



E - S(o) - £ 



?("-» 



E-E(o) -£ E-Ei'^) ~^' 

{E - £;(")) (ii; - i?(°) - e)a(t - 1') - i? - i^^°^ - e + an^-.) 

= E{n,^,) - E(n,^„) + E^°^n,^,) - E^^^n,^,) +E- S^") - £ + ({n^-„) 
= {E- i?(°))(n,_,) - (n,_,) [E - E^"^ -0+E- E^'^ - £. 
Therefore 

(ni-a) , 1 - (ni-a) 



(6.25) 



a(t-t') = 



E - {E(") + + iS E-Ei^)+i6' 
after introducing the arbitrary convergence factor, 5. 
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